Highly focused supersonic microjets 
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The paper describes the production of thin, focused microjets with velocities up to 850 m/s by 
the rapid vaporization of a small mass of liquid in an open liquid-filled capillary. The vaporiza- 
tion is caused by the absorption of a low-energy laser pulse. A likely explanation of the observed 
phenomenon is based on the impingement of the shock wave caused by the nearly-instantaneous 
vaporization on the free surface of the liquid. An experimental study of the dependence of the jet 
velocity on several parameters is conducted, and a semi-empirical relation for its prediction is devel- 
oped. The coherence of the jets, their high velocity and good reproducibility and controllability are 
unique features of the system described. A possible application is to the development of needle-free 
drug injection systems which are of great importance for global health care. 

Subject Areas: Fluid dynamics 



I. INTRODUCTION 

Transient liquid jets produced by flow focusing have 
been studied from various perspectives due to fundamen- 
tal interest and their widespread occurence (see e.g. the 
review article by Eggers and Villermaux pQ). Common 
instances include bubbles bursting at a liquid surface 
[2j [3], high amplitude Faraday waves [MS], jets gener- 
ated by the impulsive acceleration of a liquid surface [9], 
shaped charges, and others. 

In all these cases, the jet is produced when a liquid sur- 
face concave towards the gas is impulsively accelerated. 
The resulting "kinematic focusing" of the liquid, which 
converges towards the center of the curvature, is respon- 
sible for the peculiarly high-speed jets generated in this 
way. 

In shaped charges this convergence is caused by a deto- 
nation wave investing the free surface. The early history 
of these devices and the principle on which they oper- 
ate is described by Birkhoff et al. [10 , followed by many 
works, e.g., Curtis and Kelly [TT], Petit et al. [12] in more 
recent years. The focusing, however, can also be caused 
by geometrical effects simply due to the concave shape 
of the interface without reliance on liquid compressibility 
(see e.g. [IME]). 

In this study we describe yet another way for the gener- 
ation of such a jet. We use a microtube in which the im- 
pulsive acceleration of the interface is due to the abrupt 
vaporization of a small mass of liquid caused by the ab- 
sorption of a laser pulse. 

Bubble generation by the absorption of a focused laser 
pulse in the bulk of a liquid volume has been an es- 
tablished experimental technique for a long time (see 
e.g. [19j[20]). More recently, the same technique has been 
used to generate bubbles in the proximity of the free sur- 
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face of a liquid [2TH24] . In this case, strong jets are of- 
ten observed. For example, Thoroddsen et al. [24 , who 
focused a laser pulse into a drop resting on a glass sur- 
face, report jet maximum velocities of 250 m/s produced 
by small bubbles near the free surface, but without the 
ability to control the number, size and location of these 
bubbles. 

In the present work we also use a laser pulse to create 
a bubble near a free surface, but the geometry of the 
system and condition of the experiment are quite dif- 
ferent from those of earlier workers. Furthermore, the 
micro jet that we describe has extraordinary characteris- 
tics: it is created under controlled condition with good 
reproducibility; its maximum velocity reaches 850 m/s, 
namely a Mach number in air approaching 3; its diame- 
ter is typically one order of magnitude smaller than the 
diameter of the capillary tube; its generation requires 
moderate laser energies of the order of 10 /iJ, compara- 
ble to the typical energy of a laser pointer. 
A potential particularly attractive application of the 
present supersonic microjets is for needle-free drug injec- 
tion [25] . For this application, in which a liquid solution 
containing a drug is forced to penetrate human or ani- 
mal tissue through the skin, ultra-high velocities (greater 
than 200 m/s), fine scales (down to 30 /im), and good re- 
producibility and controllability are essential. 
The development of needle-free drug injection systems 
is of paramount importance in the global fight against 
the spread of disease as these systems help preventing 
contamination by needles. Several methods of microjet 
generation suitable for this application have been pro- 
posed ([26j[27]), but they all have several shortcomings: 
the production of each jet requires a significant amount 
of energy, the resulting shape of the ejected liquid mass 
tends to be diffuse rather than focused, the small size- 
nozzles are easily clogged and others. 
The microjets produced by the method described in this 
study are free from these problems. In particular, they 
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are faster and highly focused. Furthurmore, the avoid- 
ance of very thin nozzles makes the systems less suscep- 
tible to clogging, thus improving the reliability and con- 
trollability. 

The experimental setup is presented in § [TTJ followed by 
a description of jet formation and evolution in § |III| The 
quantitative dependencies of the jet speed on the various 
control parameters are discussed in § |IV[ culminating in 
a scaling law for the jet velocity as a function of the 
various parameters. The paper ends with summary and 
conclusions in 8 |Vl 
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II. EXPERIMENTAL SETUP AND CONTROL 
PARAMETERS 

Figure [TJa) shows a sketch of the whole experimen- 
tal setup, which is similar to that used by Sun et al. 
[28] . One end of the capillary tube is connected to a 
syringe pump (Model PHD 2000, Harvard Apparatus, 
USA) containing a water-based red dye. The other end 
is open to the air. A 532 nm, 6 ns laser pulse (100 mJ 
Nd:YAG laser, Solo PIV, New Wave, USA) is focused 
through a 10 x microscope objective to a point inside the 
capillary. The laser energy is monitored by means of an 
energy meter (Gentec-eo XLE4 or Gentec-eo QE12SP- 
S-MT-D0, Canada) placed behind the capillary. The 
energy absorbed by the liquid was calibrated by mea- 
suring the difference between the readings of the meter 
with the glass tube filled with the dyed and clear water. 
The jet formation was recorded using high-speed cameras 
with a frame rate of up to 10 6 fps (HPV-1, Shimadzu 
Corporation, Japan and SAl.l and Photron, USA). The 
minimum inter-frame and exposure times were 1 [is and 
250 ns respectively. This system enabled us to observe 
the capillary from below using a microscope (Axiovert 
40 CFL, Carl Zeiss, Germany). A long-distance micro- 
scope (Model K2, Infinity, USA) with a maximum mag- 
nification of 12 x was connected to the camera in order 
to capture the jet formation from the side and vary the 
field of view by adjusting the magnification. Illumina- 
tion for the camera was provided by a fibre lamp (ILP-1, 
Olympus, Japan) emitting light that passed through the 
filter protecting the camera. A digital delay generator 
(Model 555, Berkeley Nucleonics Corporation, USA) was 
used to synchronise the camera and the laser. Images 
were analysed with tracking software. 
The experimental parameters are indicated in figure [TJb) 
together with a side view of the capillary. The distance 
between the meniscus and the laser focus position is H, 
E is the absorbed laser energy, 6 is the initial contact 
angle of the liquid with the glass capillary with diameter 
D; l v is the distance offset of the laser focus with respect 
to the capillary axis along a vertical diameter, positive 
when the focus is beyond the tube axis; 1^ is the focus 
offset of the laser perpendicular to the tube axis and the 
diameter along which l v is measured. Borosilicate glass 
capillary tubes (Capillary Tube Supplies Limited, UK) of 
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FIG. 1. (a) A sketch of the experimental setup; the capillary 
tube is aligned perpendicular to the paper, (b) Definition of 
experimental parameters with a side view of the capillary in 
which the micro jet formation takes place. 



TABLE I. The range of control parameters for the initial con- 
tact angle of the liquid with the glass capillary the distance 
between the meniscus and the laser focus position H, the ab- 
sorbed laser energy E, the focus offset of the laser in the 
vertical plane l v , and the focus offset of the laser in the hori- 
zontal plane lh- The inner diameters of capillary tubes D are 
50, 200, and 500 [im. 








E([i3) 


H(fim) 


D (/im) 


l v (mm) 


l h (mm) 


Lower limit 


20° 


19 


200 


50 


-1 





Upper limit 


90° 


880 


2500 


500 


2.5 


0.1 



three different sizes were used in the experiments, with 
inner diameters D of 50, 200 and 500 /im and outer diam- 
eters of 80, 220 and 520 /im, respectively. The tip of the 
capillary tube was dipped in a hydrophobic solution (1H, 
1H, 2H, 2H-Perfluorooctyltrichlorosilane) to pin the con- 
tact line and permit a control the contact angle between 
20° and 90°. For uncoated tubes, the initial contact an- 
gle was found to be 25° ±3°, which is similar to the data 
reported by Sumner et al [29 for borosilicate glass. The 
range of the control parameters is summarised in table I. 
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III. JET FORMATION AND EVOLUTION 

An example of a supersonic jet visualized using the 
technique of dual-flash illumination by laser-induced flu- 
orescence [30 is shown in figure [2] The capillary is on 
the right, the jet tip is on the left, and the jet travels 
from right to left with a speed of 490 m/s. The tip of 
this jet has grown into an almost spherical mass which 
will eventually detach as a droplet. At the other end, the 
jet is thicker and slower. 

Figures [3]^a) and (b) show successive frames taken with 
the Shimadzu camera of the jet produced in a 500 pm 
tube by using laser energies of 365 pJ and 650 pJ. The 
bubble growth is larger in the second case and the mo- 
tion of the liquid column in the tube is more evident. The 
use of this larger capillary produces relatively slow jets 
with asymptotic velocities of ~21 m/s and 55 m/s, re- 
spectively, which facilitates the visualization. The open 
end of the tube is on the left; the other end is connected 
to the syringe pump as explained before. The bubble is 
generated near the lower wall of the capillary and mainly 
expands in the direction of the open end. In spite of 
this off-center position and asymmetric growth, the jet 
maintains an axisymmetric shape with a sharp tip. Its 
diameter is always much smaller than that of the cap- 
illary. In this instance, it is ~50 /im, about 10 times 
smaller than the capillary. The beginning of the insta- 
bility which eventually leads to droplet pinch off can be 
discerned. We have always found that, unless the bubble 
is generated less than one diameter away from the free 
surface (H < D), the asymmetry of its growth does not 
affect the axial symmetry of the jet. This observation 
supports the hypothesis that the cause of the thin jet is 
the sudden onset of a high pressure rather than the bub- 
ble expansion per se. 

Figure Qa) shows another typical sequence of jet evolu- 
tion. The corresponding history of the jet tip velocity is 
shown in Figure gb) and it is seen to be non-monotonic. 
Immediately after the laser pulse, the interface sets into 
motion (i). The maximum velocity V max occurs when 
a well formed jet just begins to appear as shown in im- 
age (ii), by which time the meniscus has lost the initial 
concave shape. This feature is observed irrespective of 
the jet speed, and is therefore compatible with the hy- 
pothesis of kinematic focusing due to the initial concave 
shape of the interface. The subsequent deceleration (hi) 
is due to surface tension as suggested by the fact that 
the deceleration is suppressed at higher velocities, i.e., 
higher Weber numbers. As the jet becomes thinner, the 
rate of new surface generation decreases, surface tension 
becomes less important, and the jet tip velocity reaches 
an asymptotic value Vj\ shortly thereafter a drop starts 
to form at the tip (iv) and it eventually detaches. 
It is worth mentioning that a very small fraction of the 
liquid receiving the initial impulse and focusing into the 
jet is in contact with the wall, and therefore its motion 
is not significantly affected by contact line phenomena. 
In any event, since the process is very rapid, such phe- 



nomena can only affect a thin liquid layer of the order of 
the viscous penetration length, which can be estimated 
as less than 5 /im. Most of the liquid is therefore essen- 
tially in an inertia-dominated, nearly inviscid motion and 
is not significantly affected by contact line phenomena. 



IV. PARAMETER DEPENDENCE 

In this section, we study the dependence of the jet 
speed on various control parameters. 

A. Contact angle 

Since the concavity of the meniscus is crucial for the 
jet formation, the first parameter we investigated is the 
initial radius of curvature Ro of the free surface. This 
quantity, which was determined by fitting a circle to the 
image of the free surface, is related to the contact angle 
6 of the liquid with the capillary by the formula 



which is valid because the effect of gravity in this sys- 
tem is negligible as shown by the smallness of the Bond 
number Bo = pgR^/j (where p is the liquid density, 7 
surface tension coefficient, and g the acceleration of grav- 
ity) typically of the order of 10 -3 . 

The initial radius of curvature is varied by minor adjust- 
ments of the liquid volume by means of the syringe. For 
these tests, the laser focus was on the tube axis and the 
energy E and the distance H between focus position and 
meniscus are kept constant at 460±20 pJ and 460±40 /im, 
respectively; experiments are carried on 500 pm diame- 
ter tubes. 

Figure |5] shows four sequences taken with different initial 
contact angles, with three snapshots per sequence. The 
shape of the jet is significantly influenced by the radius of 
curvature. These results show that the more curved the 
meniscus shape (smaller radius of curvature and lower 
initial contact angle), the higher the jet velocity due to 
the increased focusing. For increasing initial contact an- 
gles, the jet becomes thicker and less focused. Eventually, 
when the initial contact angle is 90 degrees or larger, the 
focusing is lost, no jet is formed and the liquid is only 
pushed out as a plug by the expanding bubble. 
These results confirm that the initial shape of the menis- 
cus plays an important role for the acceleration of the jet 
as discussed in section |III| For a quantitative analysis, 
we derive a scaling of the velocity increment due to the 
flow focusing as follows. 

We assume that, due to the initial high pressure pulse, 
the liquid acquires a velocity Vo normal to the initial free 
surface and the surface tension is negligible. After a short 
time At, the free surface is still spherical with a smaller 
radius Rq — VoAt. The new velocity Vq + AV acquired 



FIG. 2. Images of a supersonic micro jet generated in a 50 /xm capillary tube. The capillary is visible on the right side, the jet 
tip is shown on the left. The jet travels from right to left with a speed of 490 m/s. Time between images is 500 ns. 
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FIG. 3. Bubble growth and jet evolution after focusing the 
laser on the axis of a 500 /im capillary tube with (a) E — 
365 /xJ and H = 600 /xm and (b) E = 650 /xJ and H = 600 
/xm.The first image shows the tube at the moment the laser 
is fired. The subsequent images are taken 7/xs apart. (The 
movie is available as supplemental material [figure3.mov] ) 



due to this geometrical focusing can be estimated from 
the principle of mass conservation as 
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FIG. 4. (a) Jet evolution after focusing the laser in a 500 /xm 
capillary tube. The first image shows the tube at the moment 
that the laser is fired. The subsequent images are taken at 
later stages, each corresponding to a point in time shown in 
the adjacent figure. (The movie is available as supplemen- 
tal material [figure4a.mov].) (b) Velocity of the jet tip in a 
500 /xm tube as a function of time after the laser is fired. 



Neglecting a higher order term, we obtain the scale of 
the velocity increment as 

M-^. (3) 



Since the focusing effect is caused by the liquid par- 
ticles approaching the axis of the tube, it occurs on a 
characteristic length scale of the order of the mean dis- 
tance between the free surface of the liquid and the axis, 
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which is of the order of the tube diameter D. Thus the 
time scale At for the geometrical focusing can be esti- 
mated as 



At - 



D 

Vo' 



(4) 



This estimate is supported by the fact that, as described 
in § |III[ the acceleration due to the focusing occurs ap- 
proximately until the meniscus loses the concave shape. 
Thus, upon combining eqs. ([!]), (J3|, and Q, the increase 
in velocity due to flow focusing can be estimated as 



AV~ V cosO. 



(5) 



Thus, the velocity Vj resulting from the focusing may be 
expected to be given by 



V + AV = V o (l + 0cos0), 



(6) 



As shown in figure [6j this relation provides a reasonable 
fit to the data taken in the 500 jam tube by taking f3 
13 and Vo ~ 2.5 m/s. A more elaborated model will be 
described in Peters et al. [31 . The curve 1 — sin sug- 
gested by Antkowiak et al. [9 is also shown in figure 6] 
Although, unlike the previous model which estimates the 
asymptotic jet velocity, the latter model is aimed to de- 
scribe the initial velocity of the jet, its trend agrees with 
the present experimental data. 

If the initial velocity Vo is due to the reflection of a shock 
wave from the free surface, the strength Ap of the shock 
can be estimated as Ap = l/2pcVo, where c is the speed 
of sound in the liquid. With the previous value of Vo we 
find Ap ~ 20 atm. 

The initial velocity for the incompressible model can be 
estimated by assuming that it is acquired during a short 
time At by the action of an over-pressure Ap. If the laser 
focus is at a distance H from the free surface, this argu- 
ments leads us to HpV ~ ApAr [32 . With H = 460 
/im, Vo = 2.5 m/s, and At = 6 ns, this estimate gives 
Ap ~ 2000 atm, which is larger than the critical pres- 
sure of water, ^218 atm. For higher velocity cases, the 
incompressible pressure estimate increases considerably 



B. Distance H between laser focus and free surface 

Figure [7] shows the experimental results for the 
asymptotic jet tip velocity Vj, defined in figure |4j as a 
function of the distance H between the laser spot and 
the free surface (see figure [TJd). Two different capillary 
tubes with diameters of 200 jam and 500 /im, and two 
different energy levels E for each, 165 /iJ and 232 /iJ, 
and 305 /iJ and 458 /iJ, respectively, were used. The 
data shows that, for both tube diameters, Vj is inversely 
proportional to the distance H over a decade. This 
dependence is particularly well satisfied at the lower 
energy levels. 



It seems possible to rationalize this observation both 
on the basis of the incompressible flow focusing model 
and of the shock wave model. In the framework of the 
former, one may note that the mass of the liquid slug 
set into motion is approximately proportional to H, and 
therefore so is its inertia. In the framework of the shock 
wave mechanism, on the other hand, we recall that the 
pressure amplitude of a locally generated shock wave 
in a free fluid decreases inversely with the distance. In 
our conditions, the fluid is confined within the tube 
but the size of the initially vaporized liquid is at least 
one order of magnitude smaller than the tube diameter 
so that this geometrical attenuation may play a role. 
Another possibility is that there is viscous attenuation 
of the shock due to the no-slip condition at the tube 
wall (see e.g. [33J El]). Whatever the reason for the 
decreasing the shock pressure with increasing distance 
H, since the experiments show that the velocity is 
inversely proportional to the distance, this suggests that 
the velocity is proportional to the pressure. 



C. Absorbed energy E 

The next parameter of interest is the energy absorbed 
by the liquid. For this part of the study, the laser fo- 
cus is on the tube axis and its distance from the menis- 
cus was kept constant at 410 ± 40 /im, 390 ± 40 /im, and 
600 ± 40 jam for the 50 /im, 200 /im, and 500 jam diameter 
tubes, respectively. 

Figure [8la) shows the asymptotic jet velocity Vj vs. the 
absorbed energy for different tube diameters. The data 
are well fitted by a linear function with a positive in- 
tercept at Vj ; = 0. This feature embodies the existence 
of a threshold Eheat below which no jet is formed. The 
threshold values for the 50 /im, 200 /im, and 500 jam cap- 
illary tubes are approximately 20 /iJ, 100 /iJ, and 200 /iJ, 
respectively. Thus, the threshold value is an increasing 
function of the tube diameter D. It is interesting to note 
that, for the 50 jam capillary tube, the linear relation be- 
tween Vj and E — E^eat is preserved even when the jet 
speed becomes supersonic as shown in figure [8] (b). The 
absorbed energy for all cases is of the order of 10 9 J/m 3 , 
which is of the same order as the vaporization enthalpy 
of water at normal conditions. 

As noted before, in the framework of the shock wave 



model, the experimental results of § |IVB| suggest that 
the velocity is proportional to the pressure. This linear 
relation would then imply that the pressure of the shock 
wave at the meniscus is proportional to the absorbed en- 
ergy. From the incompressible viewpoint, one may con- 
clude that the initial impulse imparted to the liquid is 
proportional to the absorbed energy. 
The slopes of the data for the 50 /im, 200 /im, and 500 jam 
diameter tubes are 3.09, 0.26, and 0.12 m/(s-/iJ), respec- 
tively. Thus, for the same absorbed energy, the jet speed 
decreases with increasing tube diameter. The effects of 
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FIG. 5. Images showing the effect of varying the initial contact angle on the jet shape: (1) — 26°, (2) — 48°, (3) — 
67° and (4) — 74°. The first image shows the liquid in the capillary prior to the laser pulse, illustrating the initial contact 
angle and meniscus shape. The second image in each set is taken 24 /xs after the laser pulse, and the third image follows 25 fis 
later. (The movies are available as supplemental material (1) for figure5-l.mov, (2) ngure5-2.mov, (3) ngure5-3.mov and (4) 
figure5-4.mov.) 
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FIG. 6. Asymptotic jet velocity Vj as a function of the cosine 
of the initial contact angle. The solid line is a linear fit to the 
data. The dashed line is a 1 — sin fit to the data suggested 
by Antkowiak et al. [9]. 



the diameter D on the jet velocity are discussed further 
in § [TVDl 

In figure 8^b) no data points at higher Vj could be ac- 
quired because in the current configuration the glass tube 
breaks at the point where the laser focuses when the en- 
ergy E is too high. 



D. Diameter of microcapillary D 

In the study of the dependence of the jet velocity on 
the tube diameter, it was necessary to adjust the energy 
ranges depending on the diameter. On the basis of the 
results reported above, it has been concluded that, to 



FIG. 7. Asymptotic jet velocity Vj as a function of the 
distance H between the laser spot and the free surface for 
different capillary diameters and energies. The triangles rep- 
resent the data for the 200 //m tube at E — 232 /x J (A) and 
E = 165 J (A) and the circles show the results for the 500 /xm 
tube at E = 458 fi J (O) and E = 305 fi J (•). The solid and 
dashed lines are showing a -1 power law. Typical error bars 
are shown for a few data points. 



a good accuracy, Vj is proportional to (E — Eheat) and 
(1 + /3 cos 0) and inversely proportional to H. We express 
this dependency by writing 

Vj ^ m ( E - E ^ 1 + e™ sd \ ( 7 ) 

where the pre-factor f(D) embodies the dependence on 
the diameter. By trying a power law dependence, f(D) = 
foD~ a , the data is best fitted by a ~ -1.14, which may 
be considered close to -1. 

Especially the data for the 200 /im and 500 fim tubes 
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FIG. 8. (a) Asymptotic jet velocity Vj as a function of the 
energy absorbed by the liquid in the capillary tubes. The 
squares (□), triangles (a) and circles (O) represent the data 
for capillary tubes with 50 /xm, 200 /xm, and 500 /xm inner di- 
ameter, respectively. Each data point is the result of at least 
three measurements. The lines are linear fits to the data. The 
data for the 50 /im diameter tube are shown on an expanded 
scale in (b) . Each data point is the result of at least five mea- 
surements. The linear dependence on the energy is seen to 
hold also for supersonic speeds. 



FIG. 9. Pre-factor f(D) from equation (7| versus the cap- 
illary diameter. The circles, squares, and diamonds are con- 
stants obtained from the energy dependence, distance depen- 
dence, and initial contact angle dependence experiments, re- 
spectively. The value of the pre-factor increases with smaller 
capillaries. A best fit power law (dashed line) with power 
-1.14 is shown as well as a line with power -1 as shown in 
equation ^ (solid line) , which shows less agreement with the 
data, in particular for small D. 



the following empirical law is then derived: 
(E-E heat )(l + f3cos0) 



HD 



(8) 



Calculating Co by fitting the data in figure [9] we get 
Co = 0.0027 (Pa • s) _1 . It is interesting to note that Co 
has the the dimensions of an inverse dynamic viscosity. 
This feature appears to be accidental as viscous effects 
are not expected to play a significant role as noted before. 



E. Focus offsets, l v and 



agrees well with the line of slope -1, while that for 50 /am 
is above the line with slope -1. This could be attributed 
to the fact that different effects, such as laser focusing, 
start to play an important role at smaller tube sizes due 
to the higher curvature of the tube. According to the im- 
pulse pressure description suggested by Antkowiak et al. 
[9] and the analysis leading to eq. ([5|, the diameter of 
the tube should not play a role. A possible reason for 
this dependence is that the same energy delivered to a 
smaller diameter tube will give rise to a faster jet due to 
a smaller inertia. To strengthen this hypothesis it may 
be noted that, if the kinetic energy is constant for a given 
absorbed energy E, since the mass of the accelerated liq- 
uid is approximately proportional to D 2 , it follows that 
the velocity would decrease proportionally to D~ l . From 
the results that were presented in the previous sections, 



The sensitivity of the jet to the laser focus position 
relatively to the capillary axis was studied by displacing 
it in the horizontal (lh) and vertical directions (l v ) (see 
figure [lb). The results are shown in figure [To] 
Figure^lO^a) indicates that the jet velocity decreases 
monotonically for increasing 1^ all other conditions 
being held constant. We hypothesize that this trend is 
explained by a reduced volume of the vaporized liquid. 
To test this hypothesis, the jet velocity is compared to 
the initial laser-induced vaporized volume as estimated 
with the geometrical optics approximation described in 
Appendix A. As shown in figure [K^a), the jet velocity 
roughly scales with the calculated volume of the liquid 
vaporized by the laser pulse. 

Figure 10 'b) shows the jet speed versus the distance of 
the focus from the axis of the tube l v . The jet velocity 
increases with increasing l v up to l v « 1.5 mm, and then 
decreases for larger l v . The rise is likely caused by the 
larger area of the capillary tube surface illuminated by 
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FIG. 10. (a) Asymptotic jet velocity Vj (O) as a function of 
the horizontal focus displacement of the laser for the 200 /im 
tube. The black line shows the normalized vaporized liquid 
volume determined by geometrical optics approximation in 
both graphs (right axis), (b) Asymptotic jet velocity Vj as a 
function of the vertical focus displacement of the laser. The 
inserts show the capillary and the laser focus (triangle) and 
the directions of lh and l v as seen from the capillary opening. 



the laser, which increases as l v is increased by moving 
the laser towards the tube. As the vaporization always 
occurs at the capillary wall, a larger vapor mass will 
be created by increasing l v . However, when the laser 
is moved too close to the capillary, the energy per unit 
area of the laser beam decreases and eventually drops 
below the threshold required for bubble formation. This 
situation was also modeled by geometrical optics. As 
shown in figure [lO^b), the maximum jet speed and 
maximum bubble size are found at the same offset. Even 
though the overall trend is the same, the normalized 
values of the vaporized volume deviate from those of the 
jet velocity. This difference is not surprising since there 
is no reason to think that vaporized volume and jet 
growth are linearly related to each other. The analysis 
of these aspects is beyond the scope of this work; they 
will be studied in future investigations. 
The velocity increase with increasing l v shown in figure 
10 'b) is found also for other tube diameters as shown in 
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FIG. 11. Asymptotic jet velocity Vj as a function of the 
focus offset of the laser for different capillary diameters. The 
circles (O) represent measurements for the 500 /im tube; the 
triangles (a) refer to 200 /im tube; the squares (□) are for 
the 50 /xm. Each data point is the result of at least three 
measurements. 



figure [TT] where the jet speed versus the vertical focus 
offset l v is plotted for different diameters down to 50 /im. 
Jets with a velocity up to 850 m/s could be consistently 
produced with the smallest capillaries. Measurements 
beyond this velocity could not be obtained, as the 
gradual increase of the absorbing liquid volume with 
increasing l v eventually leads to such a violent bubble 
expansion that the capillary breaks. We also observed 
that, at the higher energies, the wall shear stresses 
exerted by the liquid pushed out of the capillary by 
the expanding bubble after the thin jet formation is so 
large as to shear off the capillary tip. It is possible that, 
with different materials, the velocity of the jet can be 
increased. When the jet velocity approaches the speed 
of sound in water, some limitations to further increases 
may arise due to the dominance of compressibility effects. 



V. SUMMARY AND CONCLUSIONS 

The dynamics of the high-speed microjet generated 
by a laser-induced rapid vaporization of water in a 
microtube has been studied. It has been shown that the 
jets so generated can reach speeds as high as 850 m/s 
with good controllability. 

The dependence of the jet velocity upon various con- 
trolling parameters has been investigated in a series of 
experiments the results of which have been summarised 
in the empirical relationship provided in equation (J8|. 
This equation shows the effect of the distance between 
the laser focus and the liquid meniscus at the mouth of 
the micro-tube, the absorbed laser energy, the liquid- 
tube initial contact angle and the tube diameter. 
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The jet velocity exhibits an inverse proportionality 
to the tube diameter and to the distance of the laser 
focus from the free surface, while it is proportional 
to the absorbed energy above a threshold value. The 
velocity is very strongly dependent on the curvature of 
the free surface, which is a function of the tube-liquid 
initial contact angle. It is also critically dependent 
on the amount of liquid vaporized and on its distance 
from the free surface, both of which can be varied by 
varying the position of the laser focus. In particular, 
it has been found that the offset of the laser focus 
with respect to the tube axis has a strong and very 
non-trivial effect. To elucidate the origin of this result 
we have used a geomterical-optics construction coupled 
with the Beer-Lambert law to determine approximately 
the size of the region where vaporization occurs and the 
absorbed energy. While still preliminary, the results of 
this analysis are in general agreement with the data. 
Further theoretical and numerical investigations of the 
jetting phenomenon will be addressed in Peters et at. 

ED- 

A priori results of the paper can be interpreted according 
to two different scenarios, one in which the shock wave 
generated by the nearly instantaneous vaporization plays 
a dominant role, and one in which the phenomenon 
is essentially incompressible. The estimate of the 
overpressure provided for the relatively low-velocity case 



in § |IV C| suggests that, according to the incompressible 
model, a much higher pressure is necessary to account 
for the observations as compared to the shock wave 
model. For higher velocity cases, the incompressible 
overpressure becomes even much larger than the critical 
pressure of water [31] . These considerations lend support 
to the shock wave interpretation of the results presented 
in this paper. 

The insights gained through this research and the ability 
to generate focused, controllable, and high velocity 
microjets open new doors for the realization, among 
others, of reliable needle-free drug delivery systems. 
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Appendix A: Estimate of the liquid energy 
absorption by geometrical optics 

A geometrical optics approximation was developed to 
analyze energy absorption of the laser beam energy ab- 
sorption by the liquid. The three main components of 
this calculation are: (1) splitting up the Gaussian laser 
beam into many different rays each one with a represen- 
tative energy; (2) following the path of each ray through 
the capillary, including (full or partial) wall-reflections; 
(3) using the Beer-Lambert law to model the local en- 
ergy absorbed by the liquid. This law describes the lo- 
cal irradiance /, and is given by / = Jq • 10 _es , with J 
the beam irradiance at the cylinder outer surface of the 
tube, e the absorption coefficient, s the arc length along 
the partial of the ray propagating in the liquid. The ab- 
sorption coefficient was measured as e = 84- 10 3 m _1 ; the 
local energy loss of a ray equals dl/ds. 
The liquid volumes are discretized into cells constitut- 
ing a grid around the central loss area. The local en- 
ergy lost by each ray in each cell was calculated from the 
Beer-Lambert law and add to the liquid in the cell. Sub- 
sequently, for each cell, the total energy absorbed over 
the duration of the light pulse was compared with an es- 
timate of the energy necessary for evaporation at room 
conditions, namely E^u = p(C p AT-\-AH vap )AV c , with p 
the liquid density, AT = 80 K the difference between the 
initial liquid temperature and the boiling temperature, 
C p = 4181 J/(kgK) the heat capacity at constant pres- 
sure, AH vap = 2.26 MJ/kg the vaporization enthalpy, 
and AV C the cell volume. If the energy absorbed by a 
cell exceeded Eb oi i, the cell volume was assumed to be 
vaporized instantaneously. Summation of the vaporized- 
cell volumes then provided a measure for the initial bub- 
ble volume and position. 
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